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Abstract 

We study the evolutionary conditions for Kelvin-Helmholtz (KH) instability in a Ha solar surge observed in NOAA 
AR 8227 on 1998 May 30. The jet with speeds in the range of 45-50 km s'*, width of 7 Mm, and electron number 
density of 3.83 x 10'*^ cm“^ is assumed to be confined in a twisted magnetic flux tube embedded in a magnetic field 
of 7 G. The temperature of the plasma flow is of the order of 10^ K while that of its environment is taken to be 
2 X 10^ K. The electron number density of surrounding magnetized plasma has a typical for the TR/lower corona 
region value of 2 x 10® cm“^. Under these conditions, the Alfven speed inside the jet is equal to 78.3 kms '. We 
model the surge as a moving magnetic flux tube for two magnetic field configurations; (i) a twisted tube surrounded 
by plasma with homogeneous background magnetic held, and (ii) a twisted tube which environment is plasma with 
also twisted magnetic held. The magnetic held twist in given region is characterized by the ratio of azimuthal to the 
axial magnetic held components evaluated at the flux tube radius. The numerical studies of appropriate dispersion 
relations of MHD modes supported by the plasma flow in both magnetic held configurations show that unstable against 
Kelvin-Helmholtz instability can only be the MHD waves with high negative mode numbers and the instability occurs 
at sub-Alfvenic critical flow velocities in the range of 25-50 kms^'. 
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1. Introduction 

Surges are phenomena in which dark dense mass are 
ejected in the solar atmosphere from chromospheric into 
coronal heights. Usually, they appear as straight or 
slig htly curved e le ctive structures , and they often re- 


slig htly curyed e ie ctiye structures , and tney otten re¬ 
cur dRoyj 1973allb : Syestka . 1976t Tandberg-Hanssen . 
197A Foukall 1990l) . At first, they were studied in Ha 
by Newton ( 1942 ) and by McMath and Mohlei ( 1948 ). 
They haye a typical size of 38 000-220 000 km, a trans- 
yerse yelocity of 30-200 kms ', and a lifetime of 10- 
20 min. The rotational or helical moti ons were, on oc¬ 
casions, also obserye d in surge actiyity (IGu et al.Lll994c 
Canfield et al. . 1996h . Surges were also obseryed as 


emission in He II 304 A sol ar images (IBohlin et al 
1975 : Georgakilas et al. . 1999ll . obtained by the slitless 
XUV spectrograph on a Skylab mission and with the 
Extreme Ultrayiolet Imaging Telescope on board the 
Solar and Heliospheric Observatory (SOHO), respec- 
tiyely. With deyeloping obseryational instruments and 
spacecrafts, chromospheric ejections were detected by 
using data from the 50 cm Swedish Vacuum Solar Tele¬ 
scope an d the Transition Reg ion and Coronal Explorer 
(TRACE, Brooks et al.L 2007h. as well as from t he Big 
Bear Solar Obseryatory dChen. .Bang, and Mai 2008) 
and with the Solar Ultrayiolet Measurements of Emit¬ 
ted Radiation (SUMER) spectroscopic obsery ations on 
board SOHO ( Chen, Innes. and Solanki . 2008 ). 
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Based on obseryational studies, it is now accepted 
that the driying mechanism of mass ejection in surges 
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is magnetic reconne c tion a t chromospheric heights. 
Kurokawa and Kawai ( 1993h found that Ha surges oc¬ 
cur at the very beginning phases of magnetic-flux emer¬ 
gence, and suggested that surges are produced by mag¬ 
netic reconnection between t he emerging flux and the 


pre-existing magnetic flux. Canfield et alJ ( 1996h re¬ 


ported that circumstances favorable to magnetic recon¬ 
nection are produced b y moving sa t ellite spots in a 
surge-productive region. llJddin et al.l (1201211 presented 
a multi-wavelength study of recurrent surges originated 
due to the photosph eric reconnections. Very recently, 
Chandra et al.l ( 2015 ) reported a multi-wavelength study 
of solar jets on 2010 Decem ber 11 using the So l ar Dy- 
namics Observatory (SDO, Dean Pesnell et al.l 2012h 
data. They found an increase in the amplitude of oscilla¬ 
tions close to their footpoints of the observed jets, which 
provides the evidence for the wave-ind uced reconnec¬ 


tion as a mechanism for jets triggering. ICanfleld et al. 


( 1996h also showed that a high-temperature X-ray jet 
and a cool untwisted surge can coexist lo cated side by 


side at the s ite (see Fig. 10c in their paper) . IShibata et al. 


(Il992h and lYokovama and Shibatal (Il995h succeeded to 
reproduce surge mass ejections from chromospheric 
heights by magnetic reconnection between an emerging 
flux and a pre-existing magn etic field. This num erical 
result was later confirmed bv iNindos et al.l (Il998h who 
studied the radio properties of 18 X-ray coronal je ts as 


observed by the Yohkoh SXT (ITsuneta et all 1199 Ih us 


ing Nobeyam a Radioheliog r aph 17 GHz data. From the 
SXT images, Nindos et al. ( 1998h computed the coro¬ 
nal plasma parameters at the location of the surge. At 
the time of maximum surge activity, they found elec¬ 
tron temperature - 2.8 x 10® K and emission mea¬ 
sure EM - 5.0 X 10^® cm“^, as well as derived con¬ 
straints on the ejecta electron number density, notably 
Ue < 6.1 X 10'® cm“^. Ha surg es and associated s oft X- 
ray loops were also studied by Schmieder et al. ( 19941) 
who performed simultaneous observations of NOAA 
AR 6850 on 1991 October 7, made with the MSDP 
spectrograph operating on the solar tower in Meudon 
and with the Yohkoh SXT. By measuring the volume 
emission measures of the two flaring loops (northern 
and so uthern ones) and the surg e region (mid-part of the 
surge). ISchmieder et al.l (119941) concluded that at 10:24- 
10:30 UT the temperature was (3^)x 10® K and the vol¬ 
ume emission measure was 10''^ cm“^. Assuming a vol¬ 
ume of (3-10) X10^^ cm^, they de rived an electron num 


ber density He = (3-6)xl0®cm iKavshap et al.l(120131) 


have observed a solar surge in NOAA AR 11271 using 
the SDO data on 2011 August 25, possibly triggered by 
chromospheric activity. They also measured the temper¬ 
ature and density distribution of the observed surge dur¬ 


ing its maximum rise and found an average temperature 
and a numbe r density of 2 x 10® K an d 4.17 x 10® cm“^, 
respectively. ISchmieder et al.l (Il996h studied the condi¬ 
tions for flares and surges in AR 2744 on 1980 Octo¬ 
ber 21 and 22 using observations from the Solar Max¬ 
imum Mission satellite and coordinated ground-based 
observations, which together covered a wide tempera¬ 
ture range from <10"' K to >10® K. In particular, the de¬ 
tected surge on October 22 had a total emission measure 
of 4.9 X 10"'^ cm“^ and duration of about 2000 s. The 
rough estimations of temperature and electron number 
density yielded Te ~ 10"'Kand«e ~ 10'®cm“^. Similar 
values for the temperature and electr on number density 
of Ha surges were obtained earlier bv l.Tain and Sorathia 
(Il987h from observations of a surge prominence in AR 
17212 on 1980 October 30 made at an interval of 5 s 
and 10 s in the Ha line center, through a Halle Alter of 
0.7 A passband in conjunction with a 15 cm aperture 
solar spar telescope. 

As seen, the electron number density and the tem¬ 
perature of solar surges can vary in rather wide lim¬ 


its, from ~10 


12 


cm 

i6 


and 10® K for cool Ha surges to 


~10® cm“^ and 10® K for high-temperature EUV surges. 
Since each surge is a jet in a well-defined magnetic 
flux tube, it is naturally to expect that the magnetohy¬ 
drodynamic waves propagating along the magnetized 
plasma flow can become unstable against the Kelvin- 
Helmholtz (KH) instability. It is well-known that KH 
instabilities occur when two fluids of different densities 
or different speeds flow by each other. In the solar atmo¬ 
sphere, which is made of a very hot and practically fully 
ionized plasma, the two flows come from an expanse 
of plasma erupting off the Sun’s surface as it passes 
by plasma that is not erupting. The difference in flow 
speeds and densities across this boundary sparks the in¬ 
stability that builds into the waves. When the instability 
reaches its nonlinear stage, vortices might form, recon¬ 
nection might be initiated and plasma structures might 
detach. Now, after the launch of SDO satellite, due to 
its high spatial and temporal resolution, the KH insta¬ 
bility of the coronal mass ejection reconnection outflow 
layer in the lower coro na, occurred on 2010 November 
3, has been imaged by iFoullon et al.l (12013h . Very re- 
cently, that KH instability observation was modeled by 
Zhelvazkov et al.l (l2014bl) . A concise but very good ex¬ 


ploration of KH instabilities in the solar atmosphere in 
view of thei r interpretation from observation s the reader 
can And in Tarovan and Ruderman ( 2012h . The aim 
of this study is to see whether MHD waves traveling 
along the surge jet can become unstable within its ve¬ 
locity range of 20-200 km s *. In the following Sect. |2] 
we will build up simplified models for the surge. Next 
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Sect. [3] deals with the derivation of the MHD wave dis¬ 
persion relations, while in Sect. |4] we will numerically 
analyze the dependence of the linear/threshold KH in¬ 
stability on relevant physical parameters of the surge 
and its environment. The last Sect. |5] summarizes our 
results. 


2. Surge models, basic parameters, and governing 
equations 


We explor e one o f the four Ha surges observed by 


Brooks et al.l (l2007h in the solar active region NOAA 


AR 8227 (N26°, E09°) on 1998 May 30 from 7:50 to 
16:50 UT. The electron number density derived from the 
emission measure analysis of the TRACE Fe IX 171 A 
images is equal to = 3.83 x 10*° cm“^ (the label ‘i’ 
stands for interior), and the jet velocity is of the or¬ 
der of 45-50 kms"*. Estimated chromospheric mag¬ 
netic field is Bfoot = 25 G and from the conserva¬ 
tion of magnetic flux between the reconnection region 
and the photos phere one find s that t he coronal magnetic 
field is 7-10 G. Brooks et al. (2007) claim that their ob- 
servations confirm the emerg ing flux regions model o f 


Kurokawa and Kawail (119931) and iShibata et al.l (119941) : 


moreover the clear evidence of spatiotemporal correla¬ 
tions between chromospheric (Ca II K, Ha) and coronal 
brightenings {TRACE Fe IX 171 A) indicates that the 
chromosphere is heated up to coronal temperatures at 
the surge footpoint by the energy release during mag¬ 
netic field reconnection, as shown in the numerical sim¬ 
ulation of Yokovama and Shibatal ( 1996 ). 

It is clear that the main body of the aforementioned 
surge is positioned at the TR/lower corona region. Thus, 
we can take the electron number density of the surge 
environment to be equal to «e = 2 x 10 ° cm“^ (the 
label ‘e’ stands for exterior), and its temperature rea¬ 
sonably can be Te = 2 x 10° K. Concerning the surge 
temperature, we suppose it to be equal to Ti = 10° K. 
Then with a background magnetic field - 1 G and 
a density contrast 77 = Pe/Pi = 0.052 (we assume 
that plasma densities in both media are homogeneous), 
we have the following characteristic sound and Alfven 
speed inside the surge and in the surrounding magne¬ 
tized plasma: Csi = 37 kms^', vai = 78.3 kms^' (more 
exactly 78.314 km s"*, which value determines the mag¬ 
netic field inside the jet to be equal to Bi - 7.04 G), and 
Cse = 196 kms^*, v'Ae = 341 kms“'. Thus, the plasma 
betas of the two media are respectively ySi = 0.27 and 
jSe = 0.28. Incompressible plasma is a good approxi¬ 
mation to study the KH instability, though the observed 
values are not favorable for the incompressibility. Nev- 



Figure 1: Equilibrium magnetic field geometries of a Ha solar surge. 


ertheless, in the following we consider the incompress¬ 
ible plasma inside and outside the surge. 

We model the surge as a vertically moving with a 
velocity vq cylindrical flux tube with radius a - Mjl 
(see Fig. [T]i, where M - 1 Mm is the surge width. 
Our frame of reference is attached to the TR/coronal 
plasma that implies that vq is the relative jet veloc¬ 
ity with respect to its environment. We must men¬ 
tion that because the density contrast, tj, is relatively 
high, in such a case, like in spicules, the occurrence 
of a KH instability, for instance of kink {m - 1) 
waves, becomes possible at generally high Alfven Mach 
numbers (the Alfven Mach number is defined as the 
ratio of jet velocity to Alfven speed inside the jet, 
Ma — I'o/i'Ai) and correspondingly at high critical flow 
velocities being far beyond the speeds a ccessible for 
surge s /spicules in the solar atmosphere (Zhelvazkov , 


20121 : iZhelvazkov and Zaaarashvilil 120121) . This cir¬ 


3^t] 


cumstance implies that the only possible way for 
emerging a KH instability in surges is the excitation 
of higher MHD harmonics that can become unsta- 
ble at sub-Alfvenic flow velocities in twisted tubes 
( Zaaarashvili et al. . 2010h . Here, we consider two pos¬ 
sible magnetic field geometries: (i) a moving twisted 
magnetic flux tube embedded in untwisted magnetic 
field Be (the left tube in Fig. [T]), and (ii) a moving 
twisted magnetic flux tube surrounded by plasma with 
twisted magnetic held lines (the right tube in Fig. [TJ. 
In our cylindrical coordinate system (r, ip, z) the mag¬ 
netic held has the following form: B - (o, B^{r), B^{r)^, 
and the flow profile inside the tube is vq = (0,0, v'o). In 
general, v'o can be a function of r, but we consider the 
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simplest homogeneous case. The unperturbed magnetic 
field B and the pressure p satisfy the pressure balance 
equation 


d 

d^ 


P + 




2p 


pr’ 


where A and B^, are constant. In the simpler case (left 
flux tube in Fig. [TJ the magnetic field outside the tube 
is homogeneous, - ( 0 , 0 , B^), and the solution to 
Eg. (|2]l inside the tube b ounded at the tube axis is (see 


( 1 ) 


Zaaarashvili et al.L 2014 ) 


where p is the magnetic permeability. 

As the unperturbed parameters depend on the r co¬ 
ordinate only, the perturbations can be Fourier analyzed 
with exp[-i(wr-m^-Lz)]. The equatio ns governing the 


Ptot{r ^ a) = aiI,„{Kir), 


(7) 


incom pressible plasma dynamics are (iGoossens et al 
1992h 


d^Ptot 

dr2 


£3 _d /rD\ 

^d7\c7/ 


dptot 

dr 


where I,„ is the modified Bessel function of order m and 
ffi is a constant. Transverse displacement can be written 
using Eq. (|5]) as 

_ai\ Kirl'mixir) 

^ir — 

r 




2mAajAiIm{Kir)l 


Ptot = 0. (2) 


( 8 ) 


where 



^ 2mB^ /OT „ , „\ 

Cl — , 

quency ua{ are given by 


(Q^-wi), 




fir- \ r J 


, x2ll/2 


C2 = 




4-Bl 


fi ar \ r I /ir^ 


Pi(G2-w2j -AA^oP-Jp] 

where the attenuation coefficient Ki and Alfven fre- 


(9) 

( 10 ) 


Ki — k. 


mA + k-B,, 


tOAi = 


^/PPi 


k B 


1 


OJA = 


^/pp ^/pp 

is the local Alfven frequency, 

kl - a> - k ■ Vo 


^ —Bip -t- k^ B^ 


(3) 


(4) 


and prime sign means a differentiation by the Bessel 
function argument. 

The solution to Eq. (|2]) outside the flux tube bounded 
at infinity is 


Ptot(r > a) = asK,„{k-j), 


( 11 ) 


is the Doppler-shifted frequency, and ptot is total (ther¬ 
mal -H magnetic) pressure perturbation. Radial displace¬ 
ment is expressed through the total pressure perturba¬ 
tion as 

D dptot , Cl 


where K„, is the modified Bessel function of order m and 
tte is a constant. Transverse displacement can be written 
as 

^ Q-e k-jK'^{k-j) 




C 3 dr C 3 ^“‘' 


(5) 


^ Pe ■ 


4e)’ 


( 12 ) 


The solution to this equation depends on the magnetic 
field and density profile. To obtain the dispersion rela¬ 
tion of MHD modes, we find the solutions to Eqs. © 
and (|5]) inside and outside the tube and merge the solu¬ 
tions through boundary conditions. 

3. Wave dispersion relations 

Before obtaining the wave dispersion relations for the 
two magnetic configurations, we first specify the twist 
of the magnetic flux tube to be a uniform one, that is, 

( 6 ) 


where, as before, the prime sign means a differentiation 
by the Bessel function argument, and the local Alfven 
frequency is 


WAe 


'PppI 


kzVAe- 


(13) 


Here, vac = Be/ yjpPn is the Alfven speed in the tube 
environment. 

The boundary conditions which merge the solutions 
inside and outside the twisted magnetic flux tube are the 
continuity of the radial component of the Lagrangian 
displacement 

^ir\r=a — ^sr\r=a 


B,^{0,Ar,Biz), 
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(14) 

























( 21 ) 


and the total pressure perturbation ( Bennett et al, , 19991) and 




II 

/ 

1 _ 


~ Ptot&\r=a‘> 

r=a 

(15) 

\ 

ppe - Cjj\^ ^ 


B- 

i<p £ 

Ptoli Sii 
pa 


where total pressure perturbations ptoti and ptote are 
given by Eqs. © and (HUl- Applying these boundary 
condition, after some algebra we finally derive the dis¬ 
persion relation of the normal MHD modes propagating 
along a twisted magnetic flux tube with axial mass flow 
Vo surrounded by plasma embedded in a homogeneous 
magnetic field 

F„,{Kia) - ImAcjAi! 


(q2 - 

— - ap- \ + A'^P„,{k,a)lnpi 

Pi ^ ^ 


(16) 


where, remember, Q = cu - ^ ■ vo is the Doppler-shifted 
wave frequency in the moving flux tube, and 


F,„{Kia) = 


/Cia/'X^ja) 


and Pmik^a) - 


kzaK',„{k,a) 

FfnikzCl) 


A derivation of Eq. (fTbl l starting from the basic equa- 
ti ons of ideal magnetohydrodynamics t he reader can see 


Zhelvazkov and Zagarashvifll 2012 . 


In the case when the outside magnetic field is also 
twisted (the right flux tube in Fig. [U), we consider that 


magn etic field, Be, has the form (iZaaarashvili et al 

2014 


Be - I 0, Bap—, Bfzz 


(t)!' 


(17) 


and the density is presented as p = pe(a/r)'^, so that the 
Alfven frequency 




mB^if + kzuBez 


y/J^a 


(18) 


is constant, which allows us to find an analytical solu¬ 
tion to the governing Eq. ©. The total pressure pertur¬ 
bation outside the tube is governed by the Bessel-type 
equation 


d^ptot ^ 5 dptot in^- 2\ r, 

+-^- \-T +'^e Ptot = 0, 


dF- 


dr 


where 


2 2 

n — m — 


Am B, 


2j}2 
Qif 


8 JIiSqip ^ Ae 


PPefl^ (w2 - {u- - 


(19) 


, (20) 


A solution to Eq. ( fT9l l bounded at infinity is 

cF 

Ptot(r > a) = a^—Ky{K^r), 


( 22 ) 


where v = V4 -i- rP and is a constant. 
Transversal displacement can be written as 


^er — U'e 


r 

- tte- 
a 


r [u- - K,,rK'y{K,,r) 

Fp^iaP- - -ABlpUp-lp 
2a[uF-u\^ 

a^p^ [aP- - - ABlpUp Ip 

ImBapCOAel 


aVe ("2 - - ABlpUp Ip J 


Ky(Ker). (23) 


By applying boundary conditions (O and (fTSl l. 
where the transversal displacements are given by 
Eqs. (fT^ and (|2^ . and total pressure perturbations ptoti 
and ptote, accordingly, by Eqs. © and (l22l i. we ob¬ 
tain the dispersion relation of the normal MHD modes 
propagating along a twisted magnetic flux tube with ax¬ 
ial mass flow vo surrounded by plasma embedded in 
twisted magnetic field 

( q ^ - F„(Kia) - ImAuxJ 
( q 2 - - AA^aPjppi 

Qv(Ker) - G 

~ L-H[a^ (w2 - J Qy{K,r) - g] ’ 
where 


(24) 


..\ _ '<e‘2Ky{Kea) ^ _2„ ,2 , ,2 

Qy(^ef) — , , L — a Psz\U ^Ae) 

Ky{Kf,a) ^ ' p 


A^ ^2(2 2\, 2mflBei^WAe 

H -—r-, G-2a (a> - M--—. 

Note that the left-hand sides of dispersion equations 
(Hbll and (l24l i are identical (this is not surprising), but 
the right-hand sides are completely different due to the 
very different magnetic field environments. 
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4. Numerical calculations and results 

The main goal of our study is to determine under 
which conditions the MHD waves propagating along 
the jet can become unstable. To conduct this investiga¬ 
tion, it is necessary to assume that the wave frequency 
w is a complex quantity, that is, w —» u) + iy, where 
y is the instability growth rate, while the longitudinal 
wavenumber k, is a real variable in the wave dispersion 
relation. Since the occurrence of the expected KH in¬ 
stability is determined primarily by the jet velocity, by 
searching for a critical or threshold value of it, we will 
gradually change its magnitude from zero to that criti¬ 
cal value and beyond. Thus, we have to solve the dis¬ 
persion relations in complex variables obtaining the real 
and imaginary parts of the wave frequency, or as is com¬ 
monly accepted, of the wave phase velocity Vph = 
as functions of at various values of the velocity shear 
between the surge and its environment, vq. 

Before starting the numerical job, we have to nor¬ 
malize all variables and to specified the input parame¬ 
ters. The wave phase velocity, Vph, and the other speeds 
are normalized to the Alfven speed inside the jet, VAi, 
which is calculated on using the axial magnetic fields 
Biz. The wavelength, A = 271 jk-, is normalized to the 
tube radius, a, that is equivalent to introducing a dimen¬ 
sionless wavenumber kM. For normalizing the Alfven 
frequency in the environment, aiAe, except the density 
contrast tj and the tube radius a, we have to additionally 
specify the ratio of the axial magnetic field components 
in both media, b = B^zIBi-. For our surge and its envi¬ 
ronment that ratio is equal to 0.994373 and because we 
consider the two plasmas as incompressible media, we 
shall take /? = 1. In the dimensionless analysis the flow 
speed, Vo, will be presented by the Alfven Mach number 
Ma = Vo/l'Ai. 

We first begin with the numerical solving Eq. (US. 
As we already said, the kink {m =1) mode might be in 
principal exited, but owing to the relatively high den¬ 
sity contrast, rj - 0.052, calculations show that this 
mode can become unstable at Alfven Mach numbers 
larger than 6.36, which means a critical jet velocity of 
498 km s ', which is inaccessible for solar surges. The 
excitations of unstable MHD waves with higher positive 
mode numbers m - 2 or m - 3 also requires very high 
critical plasma flow speeds beyond the upper limit of 
200 kms“* for surges. A distinctive decrease of the in¬ 
stability critical Alfven Mach number/jet speed one can 
achieve at the propagation of MHD waves with negative 
mode numbers. Numerical calculations of Eq. (O for 
the m - -3 MHD mode give in four instability win¬ 
dows on the k^a-axis whose position and width depends 




Figure 2: {Top panel) Growth rates of the unstable m = -3 MHD 
mode propagating on incompressible jets in four dilferent twisted in¬ 
ternal magnetic fields (with e = 0.052, 0.1, 0.2, and 0.4) dX p = 
0.052, b = 1, and corresponding critical Alfven Mach numbers. 
For krO = 1.8 the wavelength of the unstable m = -3 harmonic is 
^KH = 12.2 Mm, and the wave growth rate is yxH = 0.008 s“*. {Bot¬ 
tom panel) Marginal dispersion curves of the unstable m = -3 MHD 
mode for the critical Alfven Mach numbers as functions of the mag¬ 
netic field twist parameter e. At k^a =1.8 the critical jet velocity is 
v" = 45kms-‘. 

upon the magnetic held twist parameter e (see Eig. |3. 
The critical jet speeds for emergence of a KH instability 
accordingly are 50.1 kms^*, 48.2 kms^*, 46.2 kms^*, 
and 45 kms '. The two narrow windows correspond¬ 
ing to e = 0.025 and s = 0.1 are practically inapplica¬ 
ble to our surge-environment configuration: the wave¬ 
length, A - nM/k^a, of unstable m - —3 harmonics 
becomes comparable to the surge’s height; for instance, 
for k,a = 0.4 (the middle of the second instability win¬ 
dow) the wavelength of the unstable m - —3 harmonic 
at e = 0.1 is Tkh = 55 Mm. Actually only in the forth 
instability window (at e = 0.4) one can have real un¬ 
stable m - -3 MHD mode: for instance, at k^a =1.8 
the wavelength is Tkh = 12.2 Mm, and the correspond¬ 
ing dimensionless wave phase velocity growth rate is 
equal to 0.19245 which implies a wave growth rate 
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Figure 3: {Top panel) Growth rates of the unstable m = -2, m = -3, 
and m = -4 MHD modes propagating on incompressible twisted jets 
with £\ = 0.3 in a twisted external magnetic field with = 0.01 
at T] = 0.052, b = \, and critical Alfven Mach numbers equal to 
0.587, 0524, and 0.472. {Bottom panel) Marginal dispersion curves 
of the unstable m = -2, m = —3, and m = —4 MHD modes for the 
critical Alfven Mach numbers at the magnetic fields twist parameters 
£[ = 0.3 and = 0.01. The critical surge velocities of these modes 
are con'espondingly equal to 46 km s“^, 41 km s“*, and 37 km s“^. 

Tkh = 0.008 If we shift to the left at k^a =1.5 
(with maximal normalized wave phase velocity growth 
rate), the wavelength is /Ikh = 14.7 Mm, and the corre¬ 
sponding wave growth rate is a little bit higher, namely 
7kh = 0.0085 s '. The numerical study of the fluting- 
like harmonic m - —2 yields similar results; the only 
difference is that the fourth instability window locks at 
k~a - 1.6. Acceptable KH instability wavelengths and 
growth rates one can get again for e - 0.4. The criti¬ 
cal jet velocities for the m - —2 mode, however, are a 
little bit higher—they are in the range of 55-60 kms *. 
Note that normalized wave phase velocity on given dis¬ 
persion curve in the bottom panels of Figs. and [2 is 
equal to its label Ma. Therefore, the unstable perturba¬ 
tions are frozen in the flow and consequently they are 
vortices rather than waves. This observation is consis¬ 
tent with the KH instability in the hydrodynamics that 


Figure 4: {Top panel) Growth rates of the unstable m = -2, m = -3, 
and m = —4 MHD modes propagating on incompressible twisted jets 
with Si = 0.3 in a twisted external magnetic field with s^ = 0.01 
at 7} = 0.052, b = \, and corresponding critical Alfven Mach num¬ 
bers. For k-a = 2 the wavelength of the unstable m = -4 harmonic is 
^KH = 11 Mm, and the wave growth rate is yKH = 0.001 {Bottom 
panel) Marginal dispersion curves of the unstable m = -2, m = -3, 
and rn = -4 MHD modes for the critical Alfven Mach numbers at 
the magnetic fields twist parameters si = 0.3 and Sq = 0.01. The 
critical surge velocity of the ni = -4 hai'monic at k^a = 2 is equal 
to 25.2 kms“^ By contrast to the instability windows plotted in the 
bottom panel of Fig. [3] in this family of instability windows the nor¬ 
malized phase velocities are not constant. 


deals with unstable vortices. 

A much more interesting picture one obtains when 
study the more complicated case of a twisted magnetic 
flux tube surrounded by plasma embedded in a twisted 
background magnetic field. Our choice for the twist 
characteristics of the two magnetic fields (internal and 
external ones) are Si - 0.3 and = 0.01, correspond¬ 
ingly. Numerical solving Eq. (O for the three mode 
numbers m - -2, m - -3, and m = -4 gives for each 
mode number two instability windows; one of them 
with relatively high maximal growth rate, and a sec¬ 
ond window, next to the former, with one order lower 
maximal growth rate. These two families of instability 
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windows, for clarity, are presented in separate figures. 
Figs. |3] and |4] As seen from Fig. |3] a real KH instabil¬ 
ity one can observe mostly for the m - -A MHD mode 
and partly for the m - -3 harmonic. The wave growth 
rates of unstable modes are of the same order like those 
illustrated in Fig.|3 i.e., few inverse milliseconds. The 
critical flow velocities for the m - -2, m - -3, and 
m - -4 modes are correspondingly equal to 46 kms ', 
41 kms^*, and 37 kms“*. The second family of insta¬ 
bility windows, shown in Fig. 01 has two distinct pe¬ 
culiarities: first, the instability windows are shifted to 
the right-hand side of the k,a-axis, i.e., the propagation 
range of unstable MHD modes is extended, and second, 
the marginal dispersion curves are not constant—the 
normalized wave phase velocities gradually decrease 
with increasing the dimensionless wavenumber. If we 
fix the normalized wavenumber to be Ka - 2, the wave¬ 
length of the unstable m - —4 mode is Tkh = 11 Mm, 
and at lm(vph/vAi) = 0.02456 the wave growth rate 
is yKH = 0.001 s ', much lower than the observable 
growth rates in the first family of instability windows. 
Note also, that the critical jet velocity for emerging a 
KH instability now is remarkably lower—its normal¬ 
ized value is 0.322 that implies = 25. 2 kms~', i.e.. 


the ha lf of the surge speed evaluated by iBrooks et al. 


(l2007h . Thus, one can conclude that high-harmonic 
MHD modes can become unstable against the KH in¬ 
stability for accessible sub-Alfvenic velocities—this is 
more pronounced in the case when both magnetic fields 
are twisted. 


5. Discussion and conclusion 

In this paper, we have studied the condition under 
which MHD modes traveling on a Ha solar surge can 
become unstable against the Kelvin-Helmholtz insta¬ 
bility. Our model for the surge is a vertically mov¬ 
ing twisted magnetic cylindrical flux tube that might be 
surrounded by plasma embedded in homogeneous mag¬ 
netic field or by magnetized plasma with twisted mag¬ 
netic field. In each case the twist of given magnetic field 
(internal or external one) is characterized by the ratio of 
azimuthal magnetic field component at the inner surface 
of the tube to its longitudinal component. Among the 
MHD modes which propagate along the moving tube 
only the high harmonics with negative mode numbers 
can become unstable at accessible critical flow veloci¬ 
ties. One can observe practically a developing KH in¬ 
stability if only the parameter e characterizing the mag¬ 
netic field twist is large enough. For instance, at the first 
magnetic field configuration (the left flux tube in Fig. [TJ 
the m - -3 harmonic becomes unstable at critical flow 


velocity of 45 kms^' and wavelength Tkh = 12.2 Mm 
with a linear growth rate yKH = 0.008 s“' as the mag¬ 
netic field twist parameter s = 0.4. We note, that ex- 
ploring the KH instability ( Zhelvazkov et al. , 2014ah in 


a high-temper a ture s olar surge, like that observed by 
Kavshap et al.l ( 2013h . for the same mode at the same 


position on the k^a-axis and the same value of e, one 
obtains a wave growth rate of 0.033 s“' being exactly 
equal to the growth rate of the imaged KH instabil- 
ity in a coronal ma ss ejecta in the lower corona by 


Foullon et al. ( 2013h . On the other hand, the aforemen¬ 
tioned yKH = 0.008 s"' is of the same order as the 
growth rate of 0.003 s ' of vortex-shaped features along 
the interface between an erupting (dimming) region and 
the surro unding corona imaged by the SDOjAlA as re¬ 
ported by Ofman & ThompsCT ( 201 ll) . All these com¬ 
parisons allows us to believe that the KH instability 
might be imaged in surges, too. The second magnetic 
field configuration (the right flux tube in Fig. [T]i, re¬ 
veals some new aspects of the KH instability, notably 
for a fixed pair of magnetic fields twist parameters, in 
our case equal to ei = 0.3 and Sg = 0.01, for given 
high-harmonic mode one appears two instability win¬ 
dows on the kjCi-axis, next to each other. Thus, the range 
of KH instability is extended. True, in the second insta¬ 
bility windows the growth rates are much lower than in 
the standard instability windows, but nevertheless these 
weak/slow instabilities can occur. Moreover, the criti¬ 
cal flow velocity for emerging KH instability might be 
relatively lower. For example, the m - -4 harmonic 
with wavelength of 11 Mm can become unstable at jet 
speed of only =25 kms ' and its growth rate is equal 
to 0.001 s '. Solar Ha surges are generally small-scale 
eruptive events and their contributions to solar coronal 
heating due to triggered by the KH instability wave tur¬ 
bulence is modest. There are some cases, for instance as 
surges are detected in UV and EUV spectral lines, when 
the developing KH instability can bring forth a notice¬ 
able contribution to the solar corona energy budget. 
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